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ABSTRACT 

Microlensing perturbations to the magnification of gravitationally lensed quasar images are dependent on 
the angular size of the quasar. If quasar variability at visible wavelengths is caused by a change in the area 
of the accretion disk, it will affect the microlensing magnification. We derive the expected signal, assuming 
that the luminosity scales with some power of the disk area, and estimate its amplitude using simulations. We 
discuss the prospects for detecting the effect in real-world data and for using it to estimate the logarithmic 
slope of the luminosity's dependence on disk area. Such an estimate would provide a direct test of the standard 
thin accretion disk model. We tried fitting six seasons of the light curves of the lensed quasar HE 0435-1223 
including this effect as a modification to the Kochanek et al. (2006) approach to estimating time delays. We find 
a dramatic improvement in the goodness of fit and relatively plausible parameters, but a robust estimate will 
require a full numerical calculation in order to correctly model the strong correlations between the structure of 
the microlensing magnification patterns and the magnitude of the effect. We also comment briefly on the effect 
of this phenomenon for the stability of time delay estimates. 
Subject headings: gravitational lensing: strong - quasars: general 



1. INTRODUCTION 

Recent observational studies of the mi crolensing of gravi- 
tationally lensed quasars (see review by IWambsganssll2006l) 
have led to new constraints on the structure of the innermost 
regions of quasars and on the properties of foreground galax- 
ies. Because quasar accretion disks subtend angles compara- 
ble to the Einstein radii of stars in foreground galaxies, studies 
taking advantage of the dependence of microlensing magnifi- 
cations on the sourc e size have been able to estimate the sizes 
of accret i on di s ks jPoolev et al.1 [2007t lAnguita et ail 120081: 
iDai et all 120 10t iMorgan et all 1201 Oh . and multi-wavelength 
observations have ena bled estimates of the wavelength depen- 
dence of the disk size (IPoindexter et al.l2 008: Eig enbrod et alJ 
2008t iBate et al.l 120081: iFlovd et al.l 120091: iMosquera etail 
2009). Other studies have estimated the ratio of clumpy stel- 



lar mass to more smoothly distributed dark matter in lens- 
ing galaxies, providing the only constraints on stellar mass 
fractions not depe ndent on the initial mass function (e.g. 
iPooley et al.ll2009h . Time domain observations of microlens- 
ing variability are a particularly powerful tool for this type 
of investigation. For example, they have been used re- 
cently to estimate the inclinatio n of a quasar accretion disk 
(IPoindexter & Kochanekll2010al) as well as its transverse ve- 
locity relative t o the lens and the mean mass of t he stars in the 
lensing galaxy (IPoindexter & Kochanek 2010b|). 

To date, quasar microlensing has nearly always been mod- 
eled using accretion disks with time-independent sizes; this 
ought to be a good first-order approximation. But the fact that 
quasar luminosities are time-variable indicates that this model 
cannot hold in detail. Observationally, short-term variations 
in individ ual image fluxes have been observe d for a few lensed 
quasars dSchildl Il996t iKundic et al.l \l99% ISchechter et all 
120031) . leading to hypotheses involving microlensing of disks 
with fast-moving hot spots dGould & Miralda-Escudeiri997l: 
ISchechter et al.l 120031) or obscured by fast-m oving optically 
thick broad-line clouds (IWyithe & Loebll2002l) . 

In the spirit of these studies, we investigate a model of a 
time-varying accretion disk with a size proportional to some 



power of the quasar luminosity. This model introduces one 
new parameter, which is the exponent e between the disk ef- 
fective area and luminosity (i.e. A oc L c ); for simple black- 
body disks they are directly proportional (e = 1). We show that 
observations of the light curves of multiply imaged quasars 
have the potential to detect this disk size variability, and dis- 
cuss how such detections may be used to constrain the value 
of e. In Section |2] we outline the mathematical model for the 
observed light curve of a quasar image in terms of the intrin- 
sic quasar light curve. In Section[3]we study the amplitude of 
the effect for a time-variable disk by simulating realistic mi- 
crolensing scenarios. We describe in Section [4] an algorithm 
for measuring the effect in real-world light curves of lensed 
quasars, and report in Section [5] the results of applying it to 
the light curves of the quasar HE 0435-1223. In Section|6]we 
give a brief overview of our conclusions. 

2. MATHEMATICAL MODEL 

The optical and near-infrared continuum luminosity of 
bright quasars is thought to arise in a thin accretion dis k radi- 
ating as a multi-temperature blackbody (see review by iBlaesI 
120071) . If the blackbody effective temperature varies as a 
power law with radius, r e ff = Tq(R / Rq)~^ , then by integrating 
the Planck function over all radii, the total luminosity of the 
disk can be shown to be proportional to RqK(P), where K{(3) 
is a dimensionless function of the power law slope. This cal- 
culation has been used as a mea ns of estimating the disk size 
from the optical luminosity (e.g. lMorgan et al.ll2010l) . but our 
main concern is the scaling. The luminosity of a blackbody 
disk with a power-law temperature profile will scale linearly 
with its effective area Rq. 

This scaling implies that quasar variability at ultraviolet and 
optical rest wavelengths is matched by a proportional vari- 
ability in the size of the accretion disk. This is the simplest 
model for a time-variable disk size, because it does not re- 
quire the introduction of substructure in the disk or any new 
time scale. In the case of quasars that are experiencing mi- 
crolensing, the source size dependence of the microlensing 



2 



BLACKBURNE & KOCHANEK 



magnification will affect an image's light curve in phase with 
the quasar variability, boosting or suppressing it depending on 
the details of the microlens configuration. 

To quantify the effect, we start with the observed flux of a 
lensed quasar image as a function of time, 



f(t)=L(t)M+f c , 



(1) 



where L(t) is the time-dependent quasar luminosity and M 
is the effective magnification — specifically, the product of 
the macro-magnification and the microlensing magnification. 
The latter can change with time, but for this first discussion 
we will assume that it is constant on the time scales of quasar 
variability. M does, however, change with the area of the disk, 
because the microlensing magnification pattern is the convo- 
lution of the true caustic pattern with the (changing) source 
profile. The fc term denotes contaminating flux from non- 
disk sources, such as broad emission lines, the quasar host 
galaxy or the lens galaxy. The flux (and likewise the quasar 
luminosity L and disk area A) may be decomposed into a DC 
component and a small time-variable component: 



where 



f(0=(f) + Sf(t), 



(f) = (ML)(l+Q, and 
6f(t)=(M)SL(t)+(L)^SA(t) 



(2) 



(3) 



Here C=fc/ {ML) is the fractional flux contamination. Thus, 
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where e = (dlnL/dlnA)~ l is the logarithmic change in disk 
luminosity with disk area. If the luminosity of the disk is 
simply proportional to its area, then e = 1 . We have assumed 
in the second line that C <C 1, and that dlnM/dlnA is also 
small (we verify this assumption in Section |3). To convert 
this expression to magnitudes, we first separate the light curve 
into a time-averaged component (m) = -2.5 log(/) and a small 
variable component 5m(t) oc 6f / (/). Thus 



Sm(t) Srrii nt (t) ( 1 + e 
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where <5mi nt (f) oc SL/(L) is the intrinsic time-variable light 
curve of the quasar. This intrinsic variability is identical for 
all the images of a lensed quasar (once lensing time delays 
have been accounted for). Thus, for a pair of images with low 
amounts of contaminating flux, a difference in dlnM/dlnA 
will cause a difference in the amplitude of quasar variability 
seen in the two images. If the former is larger than a few per- 
cent, the latter ought to be detectable when comparing mea- 
sured light curves. Observationally, for quasar images labeled 
A and B, the intrinsic light curve measured in image B, rel- 
ative to that measured in A, would seem to be multiplied by 
the function 



C(t)=l + e 
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Figure 1. Examples of the distributions of dlnM/dlnA for a selection of 
macro-magnifications (shown at the upper right of each plot). The source 
half-light radius is fixed at 10~°' 5 = 32% of a solar mass Einstein radius, and 
the stellar mass fraction is fixed at 10%. 

3. SIMULATIONS 

We calculated the probability distribution functions of the 
value of dlnM/dlnA for several realistic microlensing pa- 
rameter values. We also created an example of a simulated 
light curve that includes the effect of a time-varying disk, in 
order to test the validity of the assumptions we made in Sec- 
tion El 

For these simulations, we used magnific ation maps created 
using the inverse ray-shooting software oflWambsgans s et"aT] 
( 119901) . The maps had a total surface mass density (or conver- 
gence) K tot roughly equal to their shear 7, divided into con- 
tributions from smooth matter and point masses. We spec- 
ify a broken power-law mass function for the microlens stars. 
Between 0.08M Q and 0.5M Q , its logarithmic slope is -1.8, 
and above 0.5 M Q it steep ens to -2.7. T his mass function 
is very similar to that of iKroupU (12001b . We cut off the 
mass function at 1.5M Q , because lens galaxies are typically 
early-type galaxies with old stellar populations. With this 
mass function, the average microlens mass is (m) = 0.247M©, 
in reasonable agreement with the recent measurement by 
iPoindexter & Kochanekl ( 1201 Obi) . 

Each map encodes the deviation in the magnification of its 
quasar image from that produced by a smooth mass distribu- 
tion as a function of the position of the source. They are 2000 
pixels on a side; this is 20 times the Einstein radius of a solar- 
mass microlens star. When projected back to the quasar, the 
side length is ~ 5 x 10 17 cm (the exact number depends on 
the redshifts of the lens and the quasar), and the pixel size is 
thus ~ 2.5 x 10 14 cm, or a few gravitational radii for a 10 9 M Q 
black hole. This is much smaller than the size of the optical 
accretion disk. 

3.1. Histograms 

For any particular quasar image, dlnM/dlnA depends on 
the microlensing caustic pattern associated with that image 
and the size of the accretion disk. One can calculate a prob- 
ability distribution for its value by first taking the difference 
of two copies of an appropriate pattern which have been con- 
volved with source profiles of slightly different areas, dividing 
the result by the difference in the logarithm of those areas, and 
finally creating a histogram of the resulting pattern. 

We created histograms for 8 macro-magnifications, 4 stel- 



EFFECT OF A TIME- VARYING ACCRETION DISK SIZE 



3 




+10 




+5 




+2 




-1.0-0.5 0.0 -1.0-0.5 0.0 -1.0-0.5 0.0 -1.0-0.5 0.0 
logr 1/2 /r Ein 

Figure 2. 50% probability bounds on |rflnM/rflnA| as a function of the 
ratio of the source half-light radius to a solar-mass Einstein radius. Macro- 
magnifications are shown in the upper right of each plot. The surface mass 
fraction in stars re*/ret t i s indicated by the symbol shape, with a legend in 
the upper left. 

lar fractions, and 3 source sizes that are typical of quasars 
lensed by foreground galaxies. We simulated finite-sized ac- 
cretion disks by convolving the magnification maps with cir- 
cular Gaussian profiles characterized by a half-light radius , 
a parameter which h as been shown (iMortonson et al.l 120051: 
iCongdon et al.ll2007h to be more important that the exact ra- 
dial profile of the source. Each difference map was created 
by subtracting a map which had been convolved with a source 
0. 1 dex smaller (in each dimension) than the chosen size from 
a map convolved with a source 0.1 dex larger. Figure Q] shows 
four such histograms, although this is only a fraction of the 96 
cases we considered. Figures[2]and[3]show the absolute value 
of d\\\M/d\\\A bounding 50% and 90% of the probability in 
each histogram, respectively. 

Clearly dlnM /din A is a small fraction of unity for most 
magnifications, stellar fractions, and source sizes. However, 
there is significant probability that it will be large enough to 
have an observable effect, particularly for small sources in 
highly-magnified saddle-point images. We stress, however, 
that even for large sources, where the first-order effects of mi- 
crolensing become very weak, there remains substantial prob- 
ability for d\nM / d\nA to have significantly nonzero value. 

3.2. Light curve simulation 

Our derivation in Section [2] assumed that the time scale for 
microlensing was long compared to the time scale for quasar 
variability. However, this is not always th e case, as quasars 
can vary on a range of time scales (e.g. iKellv et aTTl2009l : 
iKozlowski et al.l 120101) and microlens caustic crossings can 
happe n quickly for smaller source sizes (e.g. iChartas et alj 
120091) . In order to test our mathematical description, we ran 
a more detailed simulation of the microlensing of a quasar 
whose area varied linearly with its flux. 

We generated an intrinsic light curve for the source 
using the damped r a ndom walk stochastic model intro- 
duced by IKellv et ail (120091) and extensivel y conf irmed by 
IKozlowski et al.l ( 120101) and iMacLeod etail (l2010t) . An ac- 
cretion disk with a Gaussian profile and a time-averaged half- 
light radius 10 -0 ' 2 = 63% the size of a solar-mass Einstein ra- 
dius was run across two magnification maps, with a velocity 
of 1000 km s -1 (projected on the source plane). The maps 
were chosen to simulate a high-magnification, merging pair 
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Figure 3. 90% probability bounds on \d\nM/d\nA\ as a function of the 
ratio of the source half-light radius to a solar-mass Einstein radius. Macro- 
magnifications are shown in the upper right of each plot. The surface mass 
fraction in stars K*/K lot is indicated by the symbol shape, with a legend in 
the upper left. Note the change in vertical scale from Figure|2] 



of images with macro-magnifications of ±20 and stellar mass 
fractions K*/K tot = 0.1. We simulated 3000 epochs, each sep- 
arated by a day. At each time step, the area of the source was 
adjusted proportionally to the intrinsic light curve, and the 
magnification maps were convolved with the resulting source 
to find the microlensing magnifications. The simulation pro- 
duced observed light curves for the two quasar images which 
contained the imprint of the intrinsic quasar light curve as well 
as microlensing and the variation of the microlensing magni- 
fication with source area (i.e. ((t)). 

Figure |4] shows the difference between the resulting light 
curves, including the effect of a time-varying source size, and 
compares it to the same quantity with the source size held 
constant at the size associated with the mean flux. The first 
curve shows a clear echo of the intrinsic source light curve, 
modulated by £(f ) in good agreement with Equation|6] 

4. MEASURING THE EFFECT 

We adapted the tim e-delay fitting algorithm described by 
iKochanek et al.l (120061) to fit for the presence of £(f) in mea- 
sured light curves. This algorithm defines a model light curve 
consisting of the sum of a high-order polynomial s(t ) (repre- 
senting the intrinsic quasar variability) and lower-order poly- 
nomials A/i,(?) (representing microlensing fluctuations). The 
intrinsic light curve is constrained to be the same for all lensed 
images, whereas the micr olensing terms are allow ed to vary. 
We modified Equation 3 of lKochanek et al.l (120061) as follows: 



(7) 



* 2 =EE 

;=i j=i 

Here my and cry are the data and errors, respectively, and 
each Q(t) is parameterized by a low-order polynomial. (In 
practice, Q(t) is multiplied by only the terms in s(t) that are 
higher-order than N^, the order of AjU,(f), to avoid degener- 
acy.) We used a Markov Chain Monte Carlo (MCMC) tech- 
nique to minimize x 2 with respect to the time delays At, and 
the coefficients parameterizing Q(t). At each step in the chain, 
the linear parameters of s(t ) and A/Lt,-(f) were determined by 
minimizing the x 2 as in lKochanek et al.l ( 120061) . We fit the en- 
tire time series simultaneously, as opposed to breaking it into 
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Table 1 

Best-fit simulated £(f) parameters 



1000 1500 2000 
Time (days) 



3000 



Figure 4. Bottom panel: Simulated difference light curve for a pair of highly 
magnified quasar images. The dashed curve shows the difference that would 
be observed if the disk size was held fixed at its average value. The solid 
black curve includes the effects of an accretion disk with a time- varying area. 
For comparison, at the bottom, we show in gray the input quasar light curve 
scaled by a "typical" d\nM /d\nA factor of 0.1 (see Figures [2] and |5J. Top 
panel: Difference in d\aM /dlnA between the two images, as a function of 
time. This is the factor controlling the "echo" of the intrinsic light curve seen 
in the difference curve in the bottom panel. The source's average half-light 
radius is 63% of a solar-mass Einstein radius, and the stellar mass fraction is 
fixed at 10%. This simulation assumes perfect correction for the time delays 
between the images. 



seasons, because we found that in practice this gave us better 
control over "cross-talk" between the microlensing terms and 

ceo. 

We tested the efficacy of this technique by applying it to 
simulated light cur ves for 30 pairs of images. Like the sim- 
ulations in Section [3~2l these curves started with a stochastic 
intrinsic light curve, but these did not contain microlensing 
variability; we simply multiplied one of each pair of intrin- 
sic light curves by a slowly-varying curve representing £(f)- 
We deleted from each light curve 140 out of every 365 points 
to simulate six observing seasons, and fit a model to the re- 
maining data to determine the time delay and ((t ). The latter 
was modeled using Legendre polynomials to quadratic order: 

C(f) = X^=o ^'^ > '(^' wnere t is a scaled time variable running 
from -1 at the beginning of the light curve to 1 at its end. We 
used a quartic (A^ = 4) to fit the microlensing component, and 
set the parameter A to 1.0. We varied the order of the source 
polynomial A^ src between 120 and 190, and found a small de- 
creasing trend in the average x 2 P er degree of freedom for 
higher orders. However, the decrease was less than the scatter 
in the values, so we concluded that N STC = 120 was sufficient. 
On average, the algorithm successfully recovered the bj and 
time delay parameters. Table Q] compares the mean of the 30 
best-fit values to the true parameter values, and gives the scat- 
ter in the best-fit values and their typical formal uncertainty. 

These tests demonstrate that measurements of lilnM/iilnA 
in real-world lensed quasars have the potential to constrain the 
value of e, that is, the relationship between accretion disk size 
and luminosity. If the values that are measured (e.g. using 
the above method) fall in significa ntly u nlikely regions of the 
theoretical histograms (see Section UTJ, the e = 1 hypothesis, 
that is, a blackbody with power-law temperature profile, will 
be called into question. 

The presence of contaminating light from non-disk sources 
in our quasar image light curves acts as a spurious offset in 
the measured £(f) curves (see Equation |6}. It must be con- 
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trolled as much as possible in order to make robust estimates 
of the value of e. Contamination from the lens galaxy or from 
the quasar's host galaxy lensed into an Einstein ring can be 
measured and removed using high-quality ground-based or 
space-based imaging contemporaneous with some part of the 
light curve. Spectroscopic methods are necessary to remove 
contamination from non-disk regions of the quasar, such as 
emission lines or the dusty torus. However, the contamination 
from these sources will largely cancel in the absence of strong 
microlensing of these regions. 

5. APPLICATION TO HE 0435- 1223 

We have applied the adapted fitting algorithm described in 
Section [4] to a six-year /?-band light curve of the quadruply 
lensed quasar HE 0435-1223. The 332 epochs of data are the 
results of a photometric monitoring campaign between 2003 
December and 2009 April. Data at four epochs were obtained 
with SPICAM on the 3.5-meter telescope of the Apache Point 
Observatory (APO) in New Mexico. T he rest of the data w ere 
obtained with the ANDICAM camera dDePov et al.ll2003l) on 
the 1.3-meter SMARTS telescope. The first 141 epochs of 
this da ta set, including all the A PO data, were presented in de- 
tail by iKochanek et all (120061) . and the remainder of the data 
were reduced using the procedures described there. We defer 
publication of the remainder of the data, as well as a detailed 
study of the time delays and lens models, to a later study, as 
they are not our present focus. However, we show the light 
curves in Figure [5] The gray points indicate data from epochs 
that were excluded from our fit because one or more mea- 
surements differed from an initial fit to all the data enough to 
contribute more than 10 to the total \" ■ A microlensing event 
is visible in the light curve of image A between seasons 4 and 
5; because our microlensing model is a low-order polynomial 
incapable of reproducing such fine structure, the test excludes 
some points in these seasons. In hindsight, it makes sense 
to exclude these points, since our current purpose is to detect 
quasar variability, not microlensing. 

We modeled £,(f) for images A, C, and D as linear func- 
tions of time (image B was used as the reference image). For 
each (non-reference) image, the curve was parameterized by 
Q(t) = boi+bui, where i is again defined as a scaled time vari- 
able running from -1 at the beginning of the light curve to 1 
at its end. This choice allowed us to test for the presence of 
a dlnM /dlnA term even in the presence of contamination of 
the quasar light curves because the contamination C should 
not vary with time. The source light curve and microlensing 
curves were modeled as polynomials of order A^ src = 120 and 
Nfj, = 4, respectively. Including £(f ) dramatically improves the 
fit to the data from X 2 = 2036 for AW = 1 190 to X 2 = 1699 
for A'dof = 1184. We show the best-fit model light curves 
along with the data in Figure [5] and the best-fit parameters 
for the £,-(f ) curves in Table [2] The best-fit time delays are 



EFFECT OF A TIME- VARYING ACCRETION DISK SIZE 
2004 2005 2006 2007 2008 2009 




3000 



3500 



HJD- 



4000 
2450000 (days) 



4500 



5000 



Figure 5. Light curves of images A, B, C, and D (downward from top) of HE 0435-1223. For clarity, offsets of -0.35 and -0.1 magnitudes have been applied to 
images B and C, respectively. The solid lines are the best-fitting model light curves. The gray points, indicating epochs containing points that differed by more 
than y/Wa from an initial fit to all the data, were not used for the fit. 



the same within the errors whether or not C(t ) is included, and 
are bro adly consistent with those reported bv lKochanek et al] 
(120061) . in the sense that the order of arrival is preserved. In 
detail, our best-fit time delays differ from those previously 
reported by ~ 1 to 4<r (though correlations between the mea- 
sured values reduce this tension). The discrepancy is due to 
our strategy of modeling the microlensing as a single polyno- 
mial spanning all six seasons, instead of a separate polynomial 
for each season, as was done in the earlier work. When we fit 
each season separately, we recovered time delays within la 
of the previously reported values. Increasing improves the 
agreement between our delays and the reported delays, but the 
improvement is rather slow; our tests required ~ 20 before 
the discrepancy was cut in half. It is not our present goal to 
recalculate the time delays, but we discuss in Section [6] some 
potential implications of our analysis for time delay estimates. 
We performed a second fit, holding the time delays constant at 
the previously reported values, and found that the x 2 /A^j f in- 
creased to 1750/1 187, while the best- fit £(f) parameters were 
virtually unchanged. We report them in Table [2] 

The significantly non-zero linear coefficient b\ in the best- 
fit Q(t) curves for image A (see Table |2]i suggests that the 
changing accretion disk size has an appreciable effect on at 
least that image's light curve. It is worth noting that this 
image was already th ought to be affected by microlensing 
dKochanek et al.l 120061) . However, the limitations associated 
with our analysis method prevent us from placing too much 
quantitative confidence in this result. The fit requires that the 
microlensing variability, as well as the variation in £(f), be 
slow in time. It also did not give sensible results when we 
attempted to fit individual seasons because of instabilities be- 
tween the microlensing terms A/i and £(f). Finally, it takes 
no account of the correlation between the microlensing caus- 
tic pattern and its derivative with respect to source size — the 
microlensing terms and the ((t) terms were fit independently. 
Because of these limitations, and because we have analyzed 
only a single lensed quasar without carefully addressing the 
problem of contaminating light, we cannot yet comment on 
the value of e in detail. 



Table 2 

HE 0435-1223 best-fit £(f ) parameters 

Fixed delays Free delays 
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Note. — Errors represent formal 95% (Ax 2 = 4) con- 
fidence intervals. Because image B was the reference 
image, its parameters were fixed at the nominal values. 



The Bayesian Monte Carlo analysis method described by 
kl (120041) and updated t o include stellar motions by 
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Poindexter & Kochanekl (1201 Obi) addresses most of these lim- 
itations. This method simulates microlensing light curves us- 
ing the same technique as described in Section [3] but on a 
much larger scale. Simulated light curves that are found to 
be consistent with the measured data contribute to a Bayesian 
determination of parameters such as the source size and av- 
erage microlens mass. The e parameter may be added to the 
analysis as an additional parameter, and the source size al- 
lowed to vary with time when the light curves are simulated. 
This treatment allows for arbitrary microlensing variability, 
and would treat the d\nM / d\nA terms self-consistently. We 
are optimistic that this improved analysis, applied to a number 
of light curves, will allow constraints to be placed on e. 

6. CONCLUSIONS 

We have described a new effect arising from the microlens- 
ing of quasar accretion disks whose sizes vary in time. As- 
suming the change in size is proportional to some power of 
the change in luminosity of the quasar, we have derived what 
the effect would be on the light curves of lensed quasar im- 
ages (see, e.g., Equation^. The effect arises from the depen- 
dence of the microlensing magnification on the area of the 
source. This logarithmic derivative, <ilnM/<ilnA, depends 
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on the macro-magnification of the image, the surface den- 
sity of microlens stars, and the size of the source, along with 
the detailed arrangement of the microlensing caustic pattern. 
We have constructed probability distributions for the value 
of d\nM I din A for various parameter combinations, and find 
that its value can be significantly nonzero. In particular, we 
note that it does not decrease in magnitude with increasing 
source size as strongly as the microlensing magnification it- 
self does. We have also checked our mathematical model 
by constructing simulated light curves which demonstrate the 
microlensing of a time-varying disk. 

We have demonstrated how fits may be made to such 
light curves to recover £(f), defined as the difference in 
dlnM/dlnA between pairs of images. Some systematic er- 
rors are introduced by the presence of contaminating flux in 
the light curves, and we have discussed how these may be 
controlled. Applying the fitting method to the light curves of 
the four images of the lensed quasar HE 0435-1223, we find 
evidence for significant variation in this difference for image 
A relative to image B, indicating that the effect is likely oc- 
curring for this image. However, limitations of this fitting 
method, in particular the inconsistent treatment of the corre- 
lation between the microlensing and its size derivative, pre- 
vent us from saying anything quantitative. We propose, but 
have not implemented, an improve d treatmen t based on the 
Bayesian Monte Carlo method of iKochanekl ( 120041) . This 
method is fully self-consistent, and ought to be more robust. 
The one limitation of this approach is that an a priori model of 
the source variability will have to be used to control the source 
size rather than the best fit estimate of the source light curve 
generated for each trial. A determination of ((t), in particular 
for a large sample of lensed quasar images, may be used to 
test the power-law slope of the dependence of quasar lumi- 
nosity on disk size, a quantity that is predicted to be unity for 
simple blackbody disk models. 

Our estimates of the time delays of HE 0435-1223 ar e 
only roughly consistent with those of iKochanek et al.l d2006l) . 
whether or not we include the effect of C(0- Attempts to use 
the time delays of lensed quasar images for precision cos- 
mology have often been hampered by a lack of robustness in 
measured delays, though the more seri ous problem rem ains 
model degeneracies (see discussion in Kochanek 2006|). In 
our case the tension is at least partially due to our choice of a 
low-order polynomial stretching across the entire 6-year light 
curve to model the microlensing variability. Although it was 
necessary for the sake of numerical stability, this model lacks 
the freedom necessary to describe the highest-frequency mi- 
crolensing structures. This is another reason to adopt the self- 
consistent method described above. It is possible that previ- 



ously unmodeled variable source size effects are also behind 
some of the discrepancy, in which case including the £(f ) term 
would result in more accurate delay estimates. However, we 
note that we did not find significant differences in the best-fit 
time delays of HE 0435-1223 when the £(/) terms were al- 
lowed to vary versus when they were fixed at nominal values, 
nor did the best-fit £(/) terms change when the time delays 
were fixed at the previously reported values, so this remains 
speculation at best. 



This research was supported by NSF grant AST-0708082. 



REFERENCES 

Anguita, T., Faure, C, Yonehara. A.. Wambsganss, J., Kneib, J., Covone, G.. 

& Alloin, D. 2008, A&A, 481, 615 
Bate, N. F, Floyd, D. J. E., Webster, R. L., & Wyithe, J. S. B. 2008, 

MNRAS, 391, 1955 
Blaes, O. 2007, in ASP Conf. Ser. 373, The Central Engine of Active 

Galactic Nuclei, ed. L. C. Ho & J.-W. Wang, 75 
Chartas, G., Kochanek, C. S., Dai, X., Poindexter, S., & Garmire, G. 2009, 

ApJ, 693, 174 

Congdon, A. B., Keeton, C. R., & Osmer, S. J. 2007, MNRAS, 376, 263 
Dai, X., Kochanek, C. S., Chartas, G., Kozlowski, S., Morgan, C. W., 

Garmire, G., & Agol, E. 2010, ApJ, 709, 278 
DePoy, D. L., et al. 2003, in SPIE Conf. Ser. 4841, Instrument Design and 

Performance for Optical/Infrared Ground-based Telescopes, ed. M. Iye & 

A. F. M. Moorwood, 827 
Eigenbrod, A., Courbin, F., Meylan, G., Agol, E., Anguita, T., Schmidt, 

R. W., & Wambsganss, J. 2008, A&A, 490, 933 
Floyd, D. J. E., Bate, N. F, & Webster, R. L. 2009, MNRAS, 398, 233 
Gould, A. & Miralda-Escude, J. 1997, ApJ, 483, L13 
Kelly, B. C, Bechtold, J., & Siemiginowska, A. 2009, ApJ, 698, 895 
Kochanek, C. S. 2004, ApJ, 605, 58 

Kochanek, C. S. 2006, in Gravitational Lensing: Strong, Weak, and Micro, 
ed. G. Meylan, P. North, & P. Jetzer (Berlin: Springer), 91 

Kochanek, C. S., Morgan, N. D., Falco, E. E., McLeod, B. A., Winn, J. N, 
Dembicky, J., & Kelzeback, B. 2006, ApJ, 640, 47 

Kozlowski, S., et al. 2010, ApJ, 708, 927 

Kroupa, P. 2001, MNRAS, 322, 231 

Kundic, T., et al. 1997, ApJ, 482, 75 

MacLeod, C. L., et al. 2010, ApJ, 721, 1014 

Morgan, C. W., Kochanek, C. S., Morgan, N. D., & Falco, E. E. 2010, ApJ, 
712, 1129 

Mortonson, M. J., Schechter, P. L., & Wambsganss, J. 2005, ApJ, 628, 594 
Mosquera, A. M., Munoz, J. A., & Mediavilla, E. 2009, ApJ, 691, 1292 
Poindexter, S. & Kochanek, C. S. 2010a, ApJ, 712, 668 
Poindexter, S. & Kochanek, C. S. 2010b, ApJ, 712, 658 
Poindexter, S., Morgan, N., & Kochanek, C. S. 2008, ApJ, 673, 34 
Pooley, D., Blackburne, J. A., Rappaport, S., & Schechter, P. L. 2007, ApJ, 
661, 19 

Pooley, D., Rappaport, S., Blackburne, J., Schechter, P. L., Schwab, J., & 

Wambsganss, J. 2009, ApJ, 697, 1892 
Schechter, P. L., et al. 2003, ApJ, 584, 657 
Schild, R. E. 1996, ApJ, 464, 125 

Wambsganss, J. 2006, in Gravitational Lensing: Strong, Weak, and Micro, 

ed. G. Meylan, P. North, & P. Jetzer (Berlin: Springer), 453 
Wambsganss, J., Paczynski, B., & Katz, N. 1990, ApJ, 352, 407 
Wyithe, J. S. B. & Loeb, A. 2002, ApJ, 577, 615 



